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On the Positron Fraction and Cosmic-Ray Propagation Models 
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B. Burch and R. Cowsik 
Washington University in St. Louts, St. Louis, MO 63130, USA 

The positron fraction observed by PAMELA and other experiments up to ~ 100 GeV is analyzed 
in terms of models of cosmic-ray propagation. It is shown that generically we expect the positron 
fraction to reach ~ 0.6 at energies of several TeV, and its energy dependence bears an intimate but 
subtle connection with that of the boron to carbon ratio in cosmic rays. The observed positron 
fraction can be fitted in a model that assumes a significant fraction of the boron below ~ 10 GeV 
is generated through spallation of cosmic-ray nuclei in a cocoonlike region surrounding the sources, 
and the positrons of energy higher than a few GeV are almost exclusively generated through cosmic- 
ray interactions in the general interstellar medium. Such a model is consistent with the bounds on 
cosmic-ray anisotropics and other observations. 



I. INTRODUCTION 



The recent observation of the positron fraction in 
cosmic rays by PAMELA |l[ has created much ex- 
citement because of its possible connection with the 
annihilation or decay of dark matter in the Galaxy 
or with a variety of astrophysical processes (see [2j 
for references to these discussions). These suggestions 
were prompted by the recognition that the energy de- 
pendence of the positron fraction cannot be fitted by 
the comprehensive propagation model (solid line in 
Fig. []} developed by Moskalenko and Strong (M-S) 
0,1]. In Fig. [0 we show the PAMELA observations 
of the ratio, R, of the positron flux to that of the total 
electronic component in cosmic rays along with ear- 
lier observations 0-01 ■ The PAMELA measurements 
have been called anomalous as they do not conform to 
the predictions of the M-S model. Accordingly, new 
models of cosmic-ray propagation have been discussed 
(see references in Q)- 



General arguments based on cosmic-ray propaga- 
tion models indicate that the positron fraction should 
increase at high energies and asymptotically reach a 
value of ~ 0.6 at the highest energies. We note that 7- 
ray astronomy has shown that cosmic rays generated 
in the sources suffer nuclear interactions in the prox- 
imity of the sources [§4l5j . This has a strong bearing 
on the models of cosmic-ray propagation in that if 
a fraction of the B/C ratio observed in cosmic rays, 
especially at energies below ~ 10 GeV, is generated 
in a dense cocoonlike region surrounding the sources, 
then the contribution from spallation in the general 
interstellar medium would have a flat or a weak de- 
pendence on energy. Such a model is shown 
to fit the PAMELA observations and to be consistent 
with the high degree of isotropy observed in cosmic 
rays at high energies 0, fl9i - [2]] | . 



II. POSITRON FRACTION AT HIGH 
ENERGIES 

The asymptotic value of the positron fraction is es- 
timated by noting that cosmic rays observed near the 
Earth are accelerated in a set of discrete sources dis- 
tributed over the Galaxy [2^, [23[, which accelerate 
mostly electrons rather than positrons, as the Galaxy 
is made up of matter rather than antimatter. During 
the diffusive transport, the electronic component suf- 
fers loss of energy due to synchrotron radiation and 
inverse-Compton scattering on the microwave back- 
ground and other photons. As this loss increases 
quadratically with energy as bE 2 , the spectrum of the 
electronic component is sharply cut off at high en- 
ergies. Solutions to the diffusion equation [lj, [Hj|, 
which include the energy losses by electrons, yield a 
spectrum that cuts off as 



F e (E,r n ) ~ expi 



brjEE x 
4k(E x - E) 



(1) 



Here, r n is the distance to the nearest source, E x is 
the maximum energy up to which the sources acceler- 
ate electrons, the diffusion constant k as 10 28 cm 2 s _1 , 
and b as 1.6 x 10~ 3 GeV^Myr" 1 . Thus even for a 
very large value of E x , the directly accelerated elec- 
tron spectrum is cut off at Eb as 4k/ (6r 2 ) as 100 
GeV/(r n /kpc) 2 . The cutoff in the spectrum at ~ 1 
TeV observed by the HESS instrument [24[ indicates 
the presence of cosmic ray sources within ^200 pc of 
the solar system. If this is taken to be the typical 
spacing between the sources in our Galaxy, then we 
expect about 10 4 sources in this disk within a radius of 
~15 kpc [23[ ; accordingly, each of these sources need 
only to generate a very small fraction of the cosmic 
ray luminosity of the Galaxy, on the average. 

We do not expect the secondary electrons and 
positrons to exhibit such a cutoff because, unlike the 
discrete sources of primary electrons, the source func- 
tion for the secondary component extends from the 
nearest proximity to the solar neighborhood to far-off 
distances. The secondary positrons and electrons are 
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FIG. 1: The positron fraction measured by PAMELA and 
earlier measurements are shown. Gradient drifts in solar 
modulation may account for some of the difference in the 
data at E < 10 GeV The prediction of the positron 
fraction expected in the M-S model is shown as a solid line 
and in the NLB model as a dashed line. 



FIG. 2: The observed B/C ratio [H-H with the predic- 
tions from the M-S and NLB models are shown. 



III. ENERGY DEPENDENCE AT 
MODERATE ENERGIES 



generated through the 7r ± — >■ /i ± — > e ± decay chain, 
the pions being produced in high-energy interactions 
of cosmic rays with the matter in interstellar space, 
both of which arc distributed rather smoothly, with- 
out large overall gradients. Accordingly, the effects of 
the energy loss are less severe and the index of the sec- 
ondary electron spectrum at low energies is the same 
as the source spectrum, which is the same as that of 
the nucleon spectrum [25l— l27jj . At high energies, the 
secondary spectra of positrons and electrons steepens 
by one additional power. 



Fs(E) 
Fs(E) 



E -0 

£-09+1) 



for £<£ c 
for E^>E C 



(2) 



where f3 = 2.65 and E c ~ 100 - 200 GeV. 

This spectrum of the secondary electronic com- 
ponent will progressively dominate over that gener- 
ated by the discrete sources. This implies that at 
very high energies, the positron fraction simply corre- 
sponds to that in the production process in the high- 
energy collisions of cosmic rays. The fact that the 
p/n fraction in primary cosmic rays is greater than 
unity favors the production of e + over e - , reflecting 
the slightly greater production of tt + compared with 
7r~. Whereas the theoretical calculations [26|,|27| yield 
F s+ /F s _ w 1.5 — 2, the direct observations of 
produced by cosmic rays in the Earth's atmosphere 
yields F s+ /F s _ « 1.3 Then, for E > 1 TeV 



R S (E) 



F S+ (E) 



F S+ (E) + F S _(E) 



0.6. 



(3) 



We may expect that such a large value of R will be 
reached at E > 1 TeV, say beyond several TeV. 



A. Residence Time of Cosmic Rays 

There are two classes of models for cosmic-ray prop- 
agation with which to explain the measurements of the 
primary and secondary nuclei in cosmic rays as show 
in Fig. [3J In the M-S model, the secondary produc- 
tion is distributed throughout the Galaxy, and the ob- 
served decrease with energy of the ratio of secondary 
to primary nuclei is explained by an effective residence 
time of cosmic rays in the Galaxy decreasing with en- 
ergy 0, 0, ■ This decrease may be parameterized 
beyond a few GeV/n by 



t l (E) ~t L0 {E + E)- 



t L0 (E + EY 



(4) 



where Eo w 4 GeV/nucleon, and we have indicated 
A « 0.5 to reflect the full range 0.33 < A < 0.7 of the 
M-S models currently under discussion in literature. 
The value of tlo ~ 0.4 is in units of To, and E and Eq 
are expressed in GeV. Models of this class, which may 
be approximated by a leaky-box (LB) model [III [37| . 
produce a nuclear secondary to primary ratio such as 
that given by the dotted lines in Fig. [2] Note here 
that the LB model approximates the predictions of 
the M-S model also shown in Fig. [3J The second 
class of models takes explicit account of significant 
secondary production in dense regions in the vicinity 
of the primary cosmic-ray sources. Such a model may 
be realized as a nested leaky-box (NLB) [l6l. [l7|. 



B. Including Spallation in the Source Regions 

In the NLB model, it is assumed that subsequent 
to the acceleration, the cosmic rays spend some time 
in a cocoon-like region surrounding the sources, inter- 
acting with matter and generating some of the secon- 
daries, mainly at low energies. Such interactions will 
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also generate gamma rays through the ir° — > 2j de- 
cay and could be observed b y sp ace-borne gamma-ray 
telescopes like FERMI [38|, [39(1 ■ Since, according to 
the arguments summarized in Section 2, the average 
luminosity of a cosmic-ray source is rather low, their 
gamma-ray emission will be detected only in some 
favorable cases. The effective residence time in the 
cocoon, t c {E), is energy dependent, with the higher 
energy particles leaking away more rapidly from the 
cocoon. After they leak out of the cocoon into the 
interstellar medium, the cosmic rays at all energies 
up to several hundred TeV reside for an effective time 
t q before they escape from the Galaxy. In the NLB 
model, the observed energy dependence of the nu- 
clear secondary to primary ratio is fit with an energy- 
dependent leakage time t c (E) out of the cocoon and 
with a leakage time tq out of the Galaxy that is in- 
dependent, or nearly independent, of energy up to 
~1 PeV. These two contributions are depicted by the 
dashed lines in Fig. [2j and their sum is shown as a 
chain-dotted line. This shows that the residence time 
for cosmic rays inside the cocoon has a progressively 
steeper dependence on energy, and t c {E) and tq may 
conveniently be parameterized as 



r c (E) ~ TcoE*- 5l °<> E , 

T N = T C (E) + TQ. 



constant 



(5) 



Here, the lifetimes tq, tq, and tn are in Tq units 
and take on values t c q w 0.24 and tq ps 0.08 when E 
is expressed in GeV, with the parameters e = —0.01 
and 5 = 0.13. The cocoon should have a high den- 
sity so that adequate spallation might take place in 
the short amount of time that the cosmic rays spend 
around their sources. Circumstellar envelopes, dark 
clouds, molecular clouds, and giant molecular com- 
plexes are some of the candidates that may serve as 
cocoons. These have widely ran ging densities, from 
~ 10 7 cm' 3 down to ~ 10 2 cm -3 [40], and the cosmic 
rays need to spend anywhere from 10 yr to 10 5 yr in 
these regions to generate the requisite B/C ratio at 
~ 1 GeV. Since the dimensions of these regions are in- 
versely correlated with their densities, such residence 
time in the cocoon may be generated with diffusion 
constants in the range of 10 26 — 10 28 cm 2 s _1 . 

Both LB and NLB models can provide adequate 
fits to the nuclear secondary to primary ratios ob- 
served to date, even though the difference between 
them becomes progressively larger at higher energies. 
Whereas the LB models require an effective galactic 
residence time, tl(E), progressively decreasing with 
energy, the NLB models fit the data on cosmic-ray 
nuclei with a constant residence time tq at high en- 
ergies. Accordingly, LB models predict cosmic-ray 
anisotropies that increase with increasing energy, in 
conflict with the observations [4j, |l9|, |20|. In con- 
trast, NLB models predict constant anisotropies up 
to several hundred TeV, consistent with the obser- 



vations as shown in Fig. [3l To be specific, the ex- 
pected anisotropies, 8(E), are inversely proportional 
to the effective residence time of cosmic rays in the 
Galaxy. Accordingly, the anisotropy in the NLB 
model Snlb(E) is given by 



tnlb 



(E) 



tlb(E) 




(0) 



Here, 5m-s(A,E) refers to the anisotropy in the M- 
S model calculated for the two values A « 0.3 and 
A w 0.6 3, and 100 GeV refers to the energy at 
which tl(E) and tq intersect in Fig. [2] When their 
estimates are rescaled for the NLB model, according 
to Eq. 6, the expected levels of anisotro py become 
consistent with the observational limits [J, ll9M2lj . 

We can also directly estimate the anisotropy param- 
eter Snlb using the standard formula in cosmic-ray 
literature U 



"NLB 



3n\7p 
cp 



3k 
h c 



3 x 10" 



(7) 



where ho « 1 kpcw 3 x 10 21 cm is the scale height of 
the distribution of the cosmic rays above the Galac- 
tic plane and k w 10 28 cm s~ 2 is the diffusion con- 
stant of cosmic rays in the interstellar medium. This 
anisotropy is shown in Fig. [3] with an uncertainty of 
~ 200% as a gray band. Note this estimate matches 
the values of anisotropy scaled down from the M-S cal- 
culations using Eq. 6. Below ~ 1 TeV, the magnetic 
fields in the solar system, anchored at the Sun, pre- 
vent dipole anisotropies from being observed. Above 
~ 1 PeV, as we approach the knee in the cosmic-ray 
spectrum, the particles escape with increasing rapid- 
ity from the galactic volume, causing the anisotropy 
to increase. Keeping these factors in mind, we note 
that the anisotropy levels expected in the NLB model 
is consistent with the observations. 

Another difference between the two models is that 
they require different input spectra to be generated by 
the sources. To see this, let s n (E) represent the spec- 
trum of nuclei accelerated by the source by written 



s n {E) = s n0 E- a . 



(8) 



Since in the LB model these nuclei have a lifetime 
tl(E) and an interaction lifetime (cnHCTint) > the 
spectrum of cosmic-ray nuclei in the interstellar space 
becomes 



F n (E) = s n (E)- 



tl{E) 



(9) 

1 + co- vnt n H T L (E) 
« s n {E)T L {E) = s n T L0 {E + E)- A = s, M T L0 E- a - A 

for E > E . 
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FIG. 3: Measurements of the cosmic-ray anisotropy from 
various compilations 4, 19-2 J. Also plotted are the pre- 
dictions from models in Moskalenko and Strong (MS) [J] 
and the results from Eq. 6 (CB). The gray region shows 
the predicted anisotropy from Eq. 7. 

Here (cffmtnjf)" 1 is the effective mean free path for 
the loss of cosmic rays at a particular energy through 
nuclear interactions. In order to match the observed 
spectrum of cosmic rays with F n (E) = F^E^ 13 in the 
LB model, we need to set a + A = f3 rs 2.67. Thus 

s n (E) = s n0 E-^ +A ~ E~ 22 . (10) 

The calculation of the source spectrum in the NLB 
model, s nc (E) = s c oE~^, is a two-step process. The 
spectral density inside the cocoons is the product of 
s nc {E) and the leakage lifetime inside the cocoon: 

F C (E) = Snc (E)r c (E) = Snc (E)T c0 E^ 5l °s E . (11) 

These leak out into the interstellar space at a rate 
inversely proportional to the leakage lifetime from the 
cocoon so that 

F n (E) = F C {E)-^— = s nc (E) = s c0 E-t. (12) 

T C (E) 

Thus to match with the observed spectrum of cosmic 
rays with F n (E) ~ E~P we need £ = /3. This means 
that the observed cosmic rays have spectra identical 
to that accelerated by the sources, especially at high 
energies where losses due to ionization and nuclear 
interactions are small in the source regions. 

C. Derivation of the Positron Fraction 

In assessing the positron fraction in the NLB mod- 
els, we note that the secondary nuclei, such as B, 
which are generated by the spallation of primary nu- 
clei like C, have the same energy per nucleon as their 
progenitors. By contrast, positrons carry away, on the 
average, only about 5% of the energy per nucleon of 
their nuclear progenitors 0, [26|, [27j . This implies that 



even though a significant amount of B is generated 
through spallation within the cocoon, very little pro- 
duction of positrons at energies beyond 5 GeV occurs 
there. This is because the progenitors of the positrons 
with E >5 GeV should have energies beyond about 
100 GeV per nucleon and would rapidly leak out of 
the cocoon before they suffer significant nuclear in- 
teractions (see dashed line in Fig. [5]). Thus, in the 
M-S model with energy-dependent path-length distri- 
butions, and in the NLB models, we expect the source 
function for the positrons to be the same - it is simply 
proportional to the product of the observed spectrum 
of the cosmic-ray nuclei and the density of the inter- 
stellar medium and has the spectral form SqE~@. Be- 
low ~100 GeV, where the radiative energy losses are 
not significant, the observed positron fluxes would be 
the product of this source function and the residence 
time of cosmic rays in the Galaxy, tl{E) or tq, as 
relevant to the model under consideration. 

The calculation of the positron ratio in the two 
classes of models is straightforward when we note that 
its source spectrum S+a (E) in the interstellar medium 
is generated through nuclear interactions (2(| and has 
a nearly identical spectrum to that of the parent nu- 
clei, F$E~P , except that it is shifted down in energy 
by a factor n « 0.05 and multiplied by the rate of 
nuclear interactions 

S +G (E) « o^nF^rf-^E-P. (13) 

Here, erj„ is the inclusive cross section for the produc- 
tion of 7r + , which carries off a fraction n of the energy 
per nucleon of the primary cosmic-ray nucleus. The 
factor rf~ x in Eq. 13 accounts for the shift in the 
energy and the change in the energy bandwidth when 
transforming from the spectrum of the primary nuclei 
to that of the positrons. The source function is the 
same for both the M-S and NLB models. In the NLB 
model, there is an additional small contribution S+ c 
due to positron generation from nuclear interactions 
in the cocoon. Taking the expression for the spectral 
density for the nuclei in the cocoon from Eq. 11 shows 

S +C (E) « ca m n Hc s nc r t+p - 1+5log{EN) 

x fit-p-6log{E/ri) ^ Q4) 

However, this contribution is entirely negligible be- 
yond a few GeV. Therefore the steady state spectra 
F + lb(E) and F + nlb(E) ar e essentially given by the 
product of the source function and the effective life- 
time of the positrons in the Galaxy. At energies below 
^100 GeV the radiative losses are small and the effec- 
tive lifetimes in the two models are essentially given 
by the leakage lifetime tl(E) or tnlb ~ tq, respec- 
tively. Thus the positron spectra in the two models 
are given by 

F+lb = S +g {E)t l (E) 

~ E-P(E + E)- A ~E-^ +A );(15) 
F +NLB {E) - S +g (E)t g ~ E- '. (16) 
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In order to estimate the positron fraction in the 
two models, we divide the positron fluxes -F+lb and 
F+nlb by the spectral intensities of the total elec- 
tronic component in cosmic rays. A recent compila- 
tion of the observations of the total electronic com- 
ponent can be found along with a smooth fit to the 
data that includes a slight enhancement in the inten- 
sities below ~ 1 GeV, which corrects for the effects of 
modulation by the solar wind can be found in Cowsik 
and Burch [18[. It is straightforward to take ratios of 
these spectra and compare the theoretically expected 
positron fraction in the two models with the observa- 
tions shown in Fig. Q] At high energies, the spectra of 
positrons being essentially power laws in both the LB 
and NLB models, the shape of the positron fraction 
is controlled by the spectrum of the total electronic 
component. Thus, at high energies, we have 

Rn lb ~ S+{E)t g /F ± {E) 

~ £-2-65/^-3.1 _ £0.45 (17) 

Rm-s ~ S+{E)t l /F ± (E) 

~ E - 2 - 65 E- A /E- 3A ~E-°- 2 (18) 

We see in Fig. [1] that the NLB model shows the 
positron fraction increasing with energy at high en- 
ergies and the M-S model shows a declining positron 



fraction at high energies. 

IV. CONCLUSIONS 

We see that the nested leaky-box model pro- 
vides a satisfactory fit to the PAMELA observations. 
This analysis obviates the need for exotic sources 
of positrons, suggested by comparison between the 
PAMELA data and the M-S propagation model, and 
shows that the data may be accounted for by NLB 
propagation models. Since NLB models also relieve 
the anisotropy problem encountered in the LB/M-S 
class of models and qualitatively accommodate the ob- 
servations of 7T° gamma rays from regions near cosmic- 
ray sources, we conclude that the rising positron frac- 
tion observed by PAMELA is the natural result of 
cosmic-ray interactions in the interstellar medium. 
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